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The correlation between overlap intensities and level velocities has been introduced as a sensitive 
measure capable of revealing phase space localization. Previously applied to chaotic quantum sys- 
tems, here we extend the theory to near-integrable and mixed quantum systems. This measure is 
useful in the latter cases because it has the ability to highlight certain phase space structures depend- 
ing upon the perturbation used to parametrically vary the Hamiltonian. A detailed semiclassical 
theory is presented relating the correlation coefficient to the phase space weighted derivatives of the 
classical action. In the process, we confront the question of whether the Hannay-Ozorio de Almeida 
sum rules are simply extendable to mixed phase space systems. In addition, the R-scalings of the 
correlation coefficient and relevant quantities are derived for nearly integrable systems. Excellent 
agreement is found between the theory and the results for integrable billiards as well as for the 
standard map. 

PACS numbers: 05.45.Mt, 03.65.Sq 
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I. INTRODUCTION 

System response to parametric variations or pertur- 
bations is of great importance. It is a powerful exper- 
imental practice from which new information about a 
system can be extracted that is not generally available 
by other means, especially for complex systems, i.e. dis- 
ordered, interacting many-body, and/or simple chaotic 
systems. External parameters such as electromagnetic 
fields, temperatures, applied stresses, etc., are control- 
lable and often suitable for this purpose. Some specific 
examples include conductance fluctuations in quantum 
dots [l| , and variations in eigenmodes due to shape defor- 
mations of microwave cavities Q and vibrating plates Q . 
Intramolecular vibarational energy redistribution £| in 
molecules is another example wherein the technique of 
parametric variations can lead to useful insights regard- 
ing the key perturbations responsible for strong mode 
mixings. Molecular spectroscopy in external fields is an 
area of current interest and the response of a molecular 
system to external fields can be usefully analyzed from 
the parametric perspectives [5|. 

In the extreme chaotic or disordered limit, there ex- 
ist universalities in system response to perturbations [a , 
and the eigenstates respect the ergodic hypothesis 0, 
i.e. no phase space localization; there is only a scale to 
extract. On the other hand, the response of integrable or 
near-integrable systems is far richer being system depen- 
dent, and the Husimi distributions of the eigenstates are 
confined to classical tori H, . A perturbation is likely 
to give rise to subsets of eigenlevels/states that respond 
very similarly as a group, and the parametric quantities 
thereby carry far more information. For mixed systems 
both regular and chaotic trajectories form phase space 
structures and the most obvious, naive initial starting 
point would be to attempt to analyze each structure sep- 
arately. At this level of approximation, using the appro- 



priate semiclassical theory for regular and chaotic regions 
involves the respective integrable and chaotic Hannay- 
Ozorio de Almeida sum rule expressions weighted by 
the relative fractions of phase space volume. We find 
instead that the intermittent behavior of the unstable 
orbits prevents this simple picture from applying. 

A measure was introduced by one of us (S.T.) [HI that 
reveals phase space localization of quantum eigenstates. 
Originally introduced for chaotic systems, the measure 
correlates the changes in the eigenenergies due to a per- 
turbation (termed "level velocities") with the overlap 
intensities between the eigenstates and a probe state, 
which is often usefully chosen to be a coherent state. 
More recently, it has been identified as contributing to 
the scale associated with the fidelity in the weak pertur- 
bative regime [TJ. For integrable systems phase space 
localization is the norm, not the exception. There and in 
the near-integrable regime, the localization is associated 
with tori, resonance zones and stable periodic orbits. The 
correlation measure has been shown to highlight differ- 
ent features of classical phase space depending upon the 
perturbation . Several years ago Weissman and Jort- 
ner [l4| performed a similar study involving the Husimi 
distributions and parametric changes in eigenenergy lev- 
els. 

This paper is organized as follows: we begin by devel- 
oping the semiclassical theory of the correlations for near- 
integrable systems where we also remark on the more 
general theory for the level velocities in mixed systems. 
The correlations have been previously studied in highly 
excited rovibrational states in molecules where multires- 
onant Hamiltonians are applicable [T^- Section ITTT1 gives 
the semiclassical theories of the level velocities, strength 
functions and overlap intensity-level velocity correlation 
coefficients. The semiclassical theory is compared with 
various results from integrable billiards and the standard 
map which has the entire range of dynamics from inte- 
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grable to fully chaotic. We finish with a discussion of 
the similarities and distinctions between near-integrable 
and chaotic systems and some comments about mixed 
systems. 



so that the classical dynamics are essentially the same 
throughout the range, but it must also be broad enough 
to include enough eigenstates for statistical purposes. 



II. PRELIMINARIES 



III. SEMICLASSICAL DYNAMICS 



Consider a near-integrable quantum system governed 
by a smoothly parameter-dependent Hamiltonian H(X) 
with classical analog H(l, 6; A) where 



H(I,d;X)=H (I)+XH 1 (I,9) 



(1) 



Without loss of generality, the phase space volume of 
the energy surface, V , is taken to be a constant as a 
function of A. This ensures that the eigenvalues do not 
collectively drift in some direction in energy, but rather 
wander locally. The parameterized strength function is 
given by 
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S a (E, A) = — — / <■- m < 
Znfi 



iE */»/„l„-<ffW*/fi| a )^ 
Y.PcnWSi.E - E n (X)) (2) 



where p an (X) = \(a\E n (X))\ 2 . S a (E,X) is the Fourier 
transform of the autocorrelation function of a normalized 
initial state \a). Ahead S a (E, A) will denote the smooth 
part resulting from the Fourier transform of just the ex- 
tremely rapid initial decay due to the shortest time scale 
of the dynamics (zero-length trajectories). We will take 
| a) to be a Gaussian wave packet because of its ability 
to probe 'quantum phase space', but other choices may 
be useful depending on the circumstances. In one spec- 
troscopic application, we found it to be advantageous to 
take | a) as a momentum state [lflj . 

The overlap intensity-level velocity correlation coeffi- 
cient is defined as 



C a (X) = 



8E„(\) \ 
P<* n d\ I E 

a a aE 



(3) 



where er 2 and a\ are the local variances of intensities and 
level velocities, respectively. The brackets denote a local 
energy average in the neighborhood of E. It weights most 
the level velocities whose associated eigenstates possibly 
share common localization characteristics and measures 
the tendency of these levels to move in a common direc- 
tion. In this expression, the phase space volume remains 
constant so that the level velocities are zero centered 
(otherwise the mean must be subtracted), and C a (X) is 
rescaled to a unitless quantity with unit variance making 
it a true correlation coefficient. The set of states included 
in the local energy averaging can be left flexible except 
for a few constraints. Only energies where S a (E,X) is 
roughly a constant can be used or some intensity unfold- 
ing must be applied. The energy range must be small 



We develop a theory based upon semiclassical dynam- 
ics for the overlap correlation coefficient of regular states 
by examining its individual components, the level ve- 
locities and intensities. The theory for the variance of 
the level velocities involves the trace formula for near- 
integrable systems developed by Ullmo, Grinberg and 
Tomsovic 15! . The near-integrable trace formula is used 
rather than the simpler integrable trace formula of Berry 
and Tabor [0| since it correctly describes the contribu- 
tions of short, unstable periodic orbits that occur after 
a system breaks integrability. The level velocities are 
derivatives and thus far more sensitive to the perturba- 
tion than the intensities. Hence, it is more important 
to accurately describe the level velocities. The intensity 
variances for near-integrable systems are given by the in- 
tegrable result as long as the coherent state is not near a 
resonance zone. 

It turns out that a general result for the level veloc- 
ity variance in mixed systems cannot be obtained simply 
by applying the appropriate Hannay-Ozorio de Almeida 
sum rule for the effectively regular and chaotic trajec- 
tories separately. The "effective" trajectories are not 
easily defined since they do not merely correspond to 
the trajectories within the KAM islands and the chaotic 
sea, respectively. Although, it is observed that the near- 
integrable result for the level velocities is in good agree- 
ment with the calculated results before the break up of 
the last KAM torus. 



Level velocities 



B y a dapting a method employed by Berry and Keat- 
ing [13 for classical chaotic systems with the topology 
of a ring threaded by quantum flux and later generalized 
for all continuous time chaotic systems 0, 0] , the vari- 
ance of the level velocities can be obtained. The details 
for near-integrable systems are in appendix [S] where the 
variance of the level velocities is found to be 



(s\) + a 2 J 2 (s A )] 



\ dX 

x exp 



-{Max)- J 2 (s A )) 2 + Ji 2 (sa) 



2eT M V 
h Jit 



(4) 
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where the e — > limit is understood, and the amplitude 
factor is determined by the integrable system, H (I) 



Am = 



1 



2vr|M|(< i + 1 )/ 2 v ^Mj 



(5) 



The prime on the sum excludes the M = term. M 
labels the rational tori and is a <i-dimensional vector 
with positive integer components whose classical paths 
are those which at time t have returned to the same point 
on their torus after making Mi circuits of coordinate 6\ , 
M 2 circuits of coordinate 9 2 , etc. [l^. K-^i is the scalar 
curvature matrix of the energy contour and Tm is the 
period of the unperturbed orbit on the resonant torus. 
Jq(z) and J\{z) are the standard Bessel functions. By 
the Poincare-Birkhoff theorem only two periodic orbits 
survive the break-up of a rational torus. One orbit is 
stable and the other is unstable with actions S s and S u 
and stability matrices M s and M u , respectively. Hence, 
we define 



Sm (A) 



5« + S s 



and AS M (A) 



Su S s 



with s\ = AS , m(A)/?i and 



5(A) = 



K- 1 



where 



dct(A/ u - 1) 
" det(M s - 1) 



1/2 



(6) 



(7) 



(8) 



In the case k — > 1 then 5(A) — ► and the spectral stair- 
case reduces to Ozorio de Almeida's result |20j . 

The summation over rational tori can be accomplished 
by applying the Hannay-Ozorio de Almeida sum rule for 
integrable systems 



Thus, 



2e 



E 

M 



T 2 



'A 2 



V 
2^ 



dT 

J^2 



(9) 



'5S M X 1 



dX 



[J 2 ( SA ) + a 2 J 2 ( SA )] 



/ dAS M x 2 
-2eT\ 



x exp 



4 

dT 



(Jo(sa)- J 2 (sa)) 2 + J?(sa) 




For rational tori, the change of the classical action is 
proportional to the period of motion. We can think of 
this as a "ballistic" action change in contrast to the dif- 
fusive action changes of chaotic trajectories. The vari- 
ance is therefore proportional to the square of the pe- 
riod. According to the numerical calculations discussed 



ahead, ballistic action changes hold equally well for inte- 
grable and near-integrable systems, but this is not true 
for mixed phase space systems. Thus, excluding mixed 
systems, 



[J 2 ( Sa ) + 5 2 J 2 ( Sa )] 





(Jo(sx)-J 2 (sx)) 2 + Jl(sx) 



It was shown in 0] that the expression involving Bessel 
functions collapsed to the usual Berry- Tabor expression 
for most tori (with a minor adjustment for the slight dis- 
tortion of the torus, i.e. s\ — 0, dASw/dX — 0). Only 
very slowly as h — > 0, did tori other than those possessing 
the lowest order M require the Bessel function weight- 
ings. For these higher order tori, Eq. is understood 
as being just the mean square of dS-M/dX. 

Performing the integral, the variance of the level ve- 
locities becomes simply 



a% w ((E, A; ft) 



(12) 



which is only very weakly dependent on h for those low 
order tori requiring the Bessel function weightings. Thus, 
in the integrable limiting case, the ^-dependence disap- 
pears altogether. 

For mixed phase space systems, the question is whether 
a generalization of Eqs. I|9I11|I is possible; for the fully 
chaotic case, the expressions are modified by replacing T 2 
by T and the coefficients change. Is it possible to have 
fractionally weighted parts of both expressions? First, 
instead of summing over rational tori M, the sum is con- 
sidered over all periodic orbits. Given the diffusive action 
change behavior for fully chaotic systems and ballistic 
for regular systems, conceptually, the full set of action 
changes might decompose into a linear and a quadratic 
component. Continuing with this logic, Eq. @ would ap- 
ply for the periodic orbits with ballistic action changes 
and be weighted by the appropriate fraction of phase 
space volume. It is important to note a couple of sub- 
tleties here. Both a stable and unstable orbit combined to 
define {Sm, ASm}- Thus, the quadratic and linear com- 
ponents cannot be as simple as treating the stable orbits 
ballistically and the unstable orbits diffusively. Perhaps, 
at least as many orbits with positive Lyapunov exponent 
as stable orbits should be treated ballistically. Secondly, 
the key ingredient, 5m, is the average of said stable and 
unstable orbit actions. Taking the mean square of the 
sum of two orbits's actions is equal to the square of the 
sum only if the action changes are nearly equal. Exclud- 
ing the lowest order tori, this does turn out to be roughly 
true. 

Similarly, the remaining unstable orbits would be 
treated with the appropriate (chaotic) versions of 
Eqs. (|9I11|) . With this separation, the proposition would 
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be that the level velocity variance is the sum of the near- 
integrable result and the previously derived chaotic re- 
sult [H, 



4 



<;(E,\;h)f + 



gK(E,X)h d ' 



"(I-/) 



(13) 



where / is the fraction of ballistically-bchaving orbits, 
and g is a symmetry factor for the system. The quanti- 
ties C(£?, A; fi) and K (E, A) would be determined by the 
ballistic and diffusive versions of Eq. fTT|l . respectively. 
The section on the kicked rotor will address this issue. 



B. Overlap intensities 

Next, the overlap intensities are investigated and a 
semiclassical derivation of their variance is obtained. For 
an integrable system the overlap intensities are best ex- 
pressed in action-angle coordinates; details are in ap- 
pendix El and the result is 



exp 



-2(Aiy 



(14) 



M 



where in the exponential argument of the wave packet, 



each component is (AI)j = Oj(I 



The variance is 



a Gaussian weighted average of the action of the wave 
packet over the rational tori. 

The Gaussian weighting can be evaluated by replacing 
the M averaging by an average of the action variables 
over the energy surface 



exp 



-2(AI)' 



M 



V 



J exp [-2(AI) 2 /?i 2 ] 



x5(E - H (T))dl (15) 
The Hamiltonian in the numerator is expanded about I a 



giving 



i-0(I) = -H_(Ia) +<--•(_ -I~) 



(16) 



where Ho(I a ) — E and _„ = dHo(I a )/dI. Only linear 
terms are necessary in the semiclassical limit as long as 
each width component is chosen such that (Tj/h — > oo 

1/2 

as h — > 0; we take <r,j oc h 1 so that uncertainty in I 
and 9 shrink similarly as h shrinks. After performing the 
integration of the weighting over the energy shell, the 
result is 



exp 



-2(AI) ; 



M 



ttH 2 



(d-l)/2 



=(17) 



which has an ?i-scaling of /j( d-1 )/ 2 . Therefore, the lead- 
ing order ?i-scaling for the intensity variance for near 
integrable systems is 

2 _ n°j 



h (3d-l)/2 



(18) 



C. Correlations 

Finally, the semiclassical theories for the overlap in- 
tensities and the level velocities outlined above can now 
be used to construct the correlation function. Beginning 
with the covariance, it can be shown that 



d E n (X) 
dX 



C0V a (X) = (p an {X) 



2i:h d e 



V 



9N(E;X)\ 



S a (2.;A)-^ (19) 



The semiclassical expressions for the strength function 
and the parametric derivative of the staircase function 
give 



Cov a (X) 



_ £2 (d+3)/2 n(Tj _l 



/2 



\ M 

. fd AS M 

+ i 



vn 

fdS 



M 



x exp 



V dX 

T . . ia(X) 
\ ,U A Jl(SA)+ 2 
(AI) 2 2er M \\ 

n 2 n )/ E 



[-o(-a) - ia(X)Ji(s\)} 

[Jo(s\) - J 2 (sa)]J 

(20) 



Again, the diagonal approximation has been made in the 
above expression. Following arguments parallel to those 
in deriving the variances leads to the expression 



Cov a (X) = 



-2U<Tj f dT 2T/th 



9Sm 



X 



dX 



T 



exp 



~(AI) 2 



(21) 



M 



for integrable systems where th = h/e is the Heisenberg 
time. 

In order to make further progress in understanding the 
correlation function, it is necessary to know the time de- 
pendence of the M-averaged expression. It is important 
to note that the Gaussian weighting factor will not decou- 
ple from the parametric action derivative at long times 
and thus, the local average of the derivatives of the action 
is not necessarily zero. The weighted action derivatives 
can be approximated by the same method as used in the 
calculation the variance of the level velocities, so for in- 
tegrable systems 



dS 



M 



dX 



exp 



-(AI) 2 



J(dS M /dX) exp[-(AI) 2 /h 2 ]S(T - T M )dM 



M 

2 /*2i 



/ 5(T-T M )dM 



T£(A) 



(22) 



Note that the above average is linear with the period. 
The ?i-scaling of £(A) is the same as for the average 
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Gaussian weighting in the previous subsection which is 
The final result of the covariance is 



Cov a (X) 



_d/2 



(23) 



Combining this with the results above for the variances 
generates the semiclassical expression for the correlation 
function 



C a (A) 



-e(A) 



CV2(£;,A)<exp[-2(AI)V?i 2 ]) 



(24) 



M 



where the expectation value is evaluated in Eq. 117|) . For 
near-integrable systems in the semiclassical limit, the 
above result implies that Cov a (X) oc 7j( 2d_1 )/ 2 whereas 
the correlation C a (X) oc ft^ 1 )/ 4 to leading order. 



IV. RESULTS 
A. Rectangular and box billiards 

First, we consider billiards that are totally separable 
in Cartesian coordinates and thus completely integrable. 
These types of billiards have been used in the study of 
spectral rigidity |2l| . In particular, we examine a rectan- 
gular billiard with sides a, b and a box billiard with sides 
a, a, b; the sidelength a is varied. The density of states is 
kept constant by not changing the area A of the rectangle 
or the volume V of the box. The action/angle variables 
are proportional to the momentum/position coordinates 
in each direction, 



Ij = ijPj/ir and 9j = nqj/ij 



(25) 



where tj is the side length of the jth side. The Hamilto- 
nians are given by 



(26) 



The classical action of the rectangle for a closed orbit 
on a given resonant torus M is 



S m = V8mE[M 2 a 2 + M 2 (A/a) 2 } 1/2 
the period is 



T M = ^2m[M 2 a 2 + M 2 (A/a) 2 ]/E 
and the derivative of the action along the orbit is 
dS m V8mE M?a 2 ~ M%(A/a) 2 



da a [M 2 a 2 + M 2 (A/ a) 2 ] 1 / 2 

Similarly, the classical action of the box is 

5 M = V8mE[(M 2 + M 2 )a 2 + M 2 (V /a 2 ) 2 } 1 / 2 



(27) 



(28) 



(29) 



(30) 
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FIG. 1: Variance of the level velocities for rectangular and 
box billiards. The area and volume are taken to be unity, 
while the sidelength a = 1.1; the energy is 600. The upper 
panel is for the rectangle and the lower panel is for the box. 
The solid line is the quantum results while the dashed line is 
the semiclassical theory. 



the period is 



T M = j2m[M 2 a 2 



M 2 a 2 + M 2 (V/a 2 ) 2 ]/E (31) 



and the derivative of the action is 



dS M _ V8^E [Ml + Ml)a 2 - 2M$(V/a 2 ) 2 
~8a~ ~ ~ [{M 2 + M 2 )a 2 + M|(V /a 2 ) 2 } 1 / 2 



(32) 



Using the above quantities and performing the integra- 
tion in Eq. (|llf) . the variance of the action derivatives 
are 2E 2 T 2 /a 2 and 16S 2 T 2 /5a 2 for the rectangle and the 
box, respectively. Hence, the respective values of a) 
are 2E 2 /a 2 and 16£' 2 /5a 2 , which approach the variances 
of the level velocities in the semiclassical limit. Figure Q] 
demonstrates how well a) approximates the level ve- 
locities especially as T% — > 0. Note that the level velocities 
are virtually independent of Ti as predicted. 

The ?i-scaling of the intensity variance is predicted to 
be fr 5 / 2 and h 4 , respectively, for the rectangular and box 
billiards. Using Eq. I|17fl with uj a j — Tr 2 I a j /m£ 2 the 
semiclassical theory and quantum intensity variances are 
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FIG. 2: Variance of the intensities for the rectangular and box 
billiards. The area and volume are taken to be unity, while 
the sidelength a — 1.1; the energy is 600. The upper panel 
is for the rectangle and the lower panel is for the box. The 
solid line is the quantum results while the dashed line is the 
semiclassical theory. 



compared in Fig. The agreement is quite good as % — > 
0. 

The last two quantities, the covariance and correlation, 
are shown in Fig. and Fig. 0] respectively, along with 
the results from the semiclassical theory. The values of 
£(E,a) are obtained from Eq. (|22|) . For the rectangle 

(M = ^!^-(Wj (33) 

and for the box 

am z ha z \ a z a z (V/a z ) z I 

The agreement is excellent in both cases. Note that we 
have used no fitting parameters in any of the figures in 
these subsection. 




FIG. 3: Covariance for the rectangular and box billiards. The 
area and volume are taken to be unity, while the sidelength 
a = 1.1; the energy is 600. The upper panel is for the rect- 
angle and the lower panel is for the box. The solid line is 
the quantum results while the dashed line is the semiclassical 
theory. 

B. Standard map 

In this subsection we compare the semiclassical theo- 
ries to the numerical results for the standard map. The 
classical standard map is defined by 

q l+ i = qi+ Pi+i mod(l) 

Pi+i = Pi - (fc/27r)sin(27r<7i) mod(l) (35) 

where the "kicking strength" k is the varied parameter. 
The map displays the entire range of classical dynamics 
from integrable at k = to fully chaotic beyond k sa 5. 

Instead of the velocities of the eigenvalues of the Hamil- 
tonian, it is the velocities of the eigenangles for the stan- 
dard map which are of interest. The matrix elements 
of the quantum map with N discrete levels can be writ- 
ten m 

(n\U\n') = -^=exp[iir(n-n'f/N] 
yiN 

( kN \ 
xexp [i— cos[2tt(™ + a)/iV] J (36) 
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case is 



FIG. 4: Correlation for the rectangular and box billiards. The 
area and volume are taken to be unity, while the sidelength 
a = 1.1; the energy is 600. The upper panel is for the rect- 
angle and the lower panel is for the box. The solid line is 
the quantum results while the dashed line is the semiclassical 
theory. 



where n, n' = 0, . . . , iV — 1; a is a phase term which we 
set equal to zero. The effective Planck constant is given 
by h = 1/N. The eigenvalues of the propagator are 



U\i>j 



exp(^—i(f)j)\ij)j) with real eigenangles <j)j. The 



"level velocities" are given by the Hellmann-Feynman 
theorem 



^ = 2^i|cos(27r«)|Vi) 



(37) 



The smoothed density of states for the eigenangles \2c 



is 



N 1 

d e (cf>) = — + — > exp(m^)TrC7" exp(-ne) (38) 

Z7T Z7T ^ — ' 



where the sum excludes n = and the trace of the prop- 
agator in the integrable regime is 



Tr U n = V J— exp(2mNS M ) 
V in 

M 



(39) 



Nd> i 



i 



oo 1 

E E 

n—~ oo IVI 




; exp(27rzA^5 f ]vi) exp(— ne) 



(40) 



The resonant tori M including repetitions for the stan- 
dard map are given by a constant momentum equal to 
the rational fraction m/n where m = 0, . . . , n — 1 and n 
is the period of the torus. The derivative of the action 
with respect to the kicking strength is 



dS 



M 



1 



dk 



cos(27rgi) 



(41) 



For the integrable case (fc = 0), the periodic points are 
Qi = qo + m- i/n mod(l) and hence, OSm./ dk = for all 
resonant tori. Thus, by Eqs. (|llf> and (|12l) the variance 
of the level velocities is zero. 

In the near-integrable regime, the variance of the eige- 
nangles is approximately given by the result 



4ir 2 (N 2 



(42) 



Note that the eigenangles are divided by an effective h 
making the variance of the eigenangles quadratic in N. 
We do not have an analytic formula for but it can be 
obtained numerically from the quadratic period depen- 
dence of the variance of the action derivatives averaged 
over the tori. For k = 0, it is straightforward to calcu- 
late analytically the positions of the stable and unstable 
periodic orbits that arise as soon as integrability is bro- 
ken. Curiously, using these points as initial conditions 
for the periodic orbits throughout the near-integrable 
regime, roughly k G [0, 2], gives results within a few per- 
cent of that found by locating the true periodic orbits. 
In order to make long time calculations, this compromise 
is indispensable and does not harm the approximation 
noticeably. In Fig. |5j the near-integrable theory is shown 
to give excellent agreement with the quantum eigenangle 
velocity variance even though there exists a significant 
proportion of unstable orbits for the larger values of k 
shown that are being treated as though they were stable. 

In the mixed regime, the corresponding proposition of 
Eq. O is 

„l*^((k)fN 2 +g^(l-f)N (43) 

where 5 = 4 and D(k) is the classical action diffusion con- 
stant. It was shown previously [24| that a direct method 
based on classical orbits gave a good correspondence be- 
tween the theory and quantum values for 



4^ 

ci(fc)JV + c 2 (fc)iV 2 



V(q )(V(q( n \qo,Po))}ndp dq 



(44) 



Integrating over the angle, the smoothed spectral stair- 



where V(q) — cos(27r<;) but no understanding was given 
of the two coefficients. Equation (|43|l gives one plausible, 
but naive, interpretation. 



8 




FIG. 5: The variance of the level velocities for the standard 
map (N = 500). The solid line is the quantum results while 
the dashed line is the semiclassical theory using the near- 
integrable results of Eq. 11421 . 



FIG. 6: Histogram of finite time Lyapunov exponents in the 
standard map in the mixed regime. The parameters are k = 
2.7 and T = 500. The values of [t, were obtained by a least 
squares fit of the natural logarithm of Eq. H45II . 



In order to test numerically the separation of diffu- 
sive and ballistic contributions in the mixed regime, it is 
impossible to work with the periodic orbits. They prolif- 
erate exponentially and must be found numerically. The 
time over which this can be done by brute force calcula- 
tion is too short for our purposes. Instead, a Monte Carlo 
technique is needed. We have devised a method that be- 
gins with a uniform sampling of initial conditions in the 
phase space, and relies on the properties of finite time 
Lyapunov exponents j2^. Although, a stable periodic 
orbit has a vanishing Lyapunov exponent, at a partial 
period of its motion, the trace, Q, of its stability matrix 
can grow as a power of the propagation time. This can 
give a stable orbit the appearance of being unstable de- 
pending on circumstances. For example, constructing a 
surface of section based on whether \Q\ < 2 fails miser- 
ably. If instead, one seeks the least square deviation of 
the time development of Q from the form 



Q = aT 1 exp fiT 



(45) 



then the estimate of the Lyapunov exponent can be used 
to judge whether an orbit is in the stable or unstable part 
of the phase space; i.e. the construction of the correct 
surface of section is recovered in this way. 

Figure E| gives a histogram of the finite time Lyapunov 
exponents, /i, calculated using Eq. I|45l) for the orbits in 
this regime. The orbits are separated into three rough 
categories depending upon their stabilities as mentioned 
in the previous section. The broadened 5-function feature 
near zero represents the stable orbits that have action 
derivatives with a quadratic dependence upon time. The 
set of most highly unstable or chaotic trajectories give 
an approximate log-normal distribution |25| . In between 
these two features lie the weakly unstable orbits. It turns 
out that the coefficient converges well if averaged 
over the stable orbits and up to a similar number of the 
weakly unstable orbits. However, it also turns out that 



when compared to the coefficient 02(h) of Eq. I)44[l. the 
fraction / deduced by their comparison is inconsistent 
with the relative fraction used to derive £(fc) for k values 
in the range [1.5,4.0]; see appendix IO 

A second inconsistency with the simplistic separation 
problem is that the remaining strongly unstable orbits do 
not show a purely diffusive behavior. They show both a 
quadratic and linear dependence in their action derivative 
variances; i.e. the numerical evaluation of D(k) is sus- 
pect. It could also be evaluated according to its asymp- 
totic expression |24j , but again requires a large value of / 
to be consistent with the coefficient c\(k) unless k exceeds 
roughly 5 or so. At that point / « and the system is 
basically in the fully chaotic regime. We attempted sev- 
eral schemes of separating out orbits according to ballis- 
tic or diffusive behaviors, none of which worked and are 
outlined in appendix [U] 

In spite of the lack of apparent separability, the sim- 
plistic expression, E q. (1441) works quite well. We recal- 
culate the results of [24| taking into account the Heisen- 
berg time associated with the mean eigenangle spacing in 
Fig- El As a function of kicking strength k the agreement 
between the semiclassical theory and the variance of the 
eigenangle velocities is quite good throughout the entire 
transition from integrability to full chaos. 

Next we turn to the overlap intensities. Since the en- 
ergy is not conserved for maps, Eq. (|17|l does not apply. 
The variance is obtained directly from Eq. I|14|l where the 
average is over all resonant tori; it becomes just the Gaus- 
sian normalization factor. The width of the Gaussian 
wave packet in Cartesian coordinates is taken to be pro- 
portional to the square root of U, a = 0yh — /3/y2nN. 
Thus, for k = the variance of the intensities is 



(46) 



This expression does not take into account the reflection, 
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FIG. 7: The variance of the level velocities for the standard 
map (N = 500). The solid line is the quantum results while 
the dashed line the classical phase space average of Eq. II44I . 
The time, n, is taken to be N/g = 125. 



FIG. 9: Covariance of the standard map for a section of phase 
space at q a — 0.25. The parameters are TV = 500 and k = 
0.01. The solid line is the quantum results while the dashed 
line is the semiclassical theory. 
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FIG. 8: Variance of the intensities of the standard map for 
a section of phase space at q a = 0.25. The parameters are 
TV = 500 and k = 0.01. The solid line is the quantum results 
while the dashed line is the semiclassical theory. 



FIG. 10: Correlation of the standard map for a section of 
phase space at q a — 0.25. The parameters are N — 500 and 
k — 0.01. The solid line is the quantum results while the 
dashed line is the semiclassical theory. 



(q,p) — > (1 — Ojl—p), and time reversal symmetries of the 
kicked rotor [26j |. Thus, we symmetrize the wave packet 
accordingly. 

Figure [S] shows the good agreement between the semi- 
classical theory and the intensity variance for the inte- 
grable regime at k = 0.01. The only exceptional cases 
are the tori at p — and p — 0.5 which happen to be self- 
symmetric and on the lowest order resonances. Whereas 
for k — 0, they would lead to double the intensity vari- 
ance and be accurately accounted for once the symmetry 
was fully incorporated, at k = 0.01 they are strongly af- 
fected by their location on low order resonances. As we 
do not have the generalizations necessary for those cases, 
we avoid the resonance zones instead. 

Finally, we turn to the correlations. The covariance 
between the intensities and level velocities is given by 



Eq. ill'.'il divided by Ti — l/2irN since the level velocities 
in the standard map are actually eigenangle velocities. 
After symmetrizing the wave packet, Fig. [5] compares the 
semiclassical prediction with the quantum results. Also, 
the correlation is compared in Fig. ^3 In both cases 
the agreement between the semiclassical theory and the 
quantum results is good. 

In the near-integrable regime, the covariance and cor- 
relation is strongest for resonance zones; see Fig. II II and 
Ref. . The quantizing tori are deforming at a greater 
rate inside the resonance while the tori outside the reso- 
nance are only slightly perturbed. Thus, the associated 
level velocities are larger for this area of phase space. 
Also, note that there is a dip in the correlation near the 
stable periodic orbit at the center of the resonance zone. 
The dip is only present in the correlation and not in the 
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Finally, the correlation between the level velocities and 
the intensities has been shown to highlight resonance 
zones in near-integrable systems. The resonance zones 
exhibit greater localization since the motion is bounded 
within the separatrix. Investigating the correlation and 
covariance separably may yield more information about 
the system. An example is the dip in the correlation 
that occurs near the center of the resonance zone due to 
the stable periodic orbit; the dip is not present in the 
covariance. 



FIG. 11: Correlation of the standard map in the near- 
integrable regime. The parameters are TV = 500 and k — 0.5. 
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APPENDIX A: LEVEL VELOCITIES 



V. DISCUSSION 

We have derived semiclassical expressions of the corre- 
lation of level velocities and overlap intensities in near- 
integrable systems. Their ?i-scaling is very different than 
previously derived in the chaotic limit |l9j | where the level 
velocity variance is 



d-l 



2 t. 

a E oc n 

and the intensity variance IS 



2 *. 

a„ <x n 



2d— 1 



(47) 



(48) 



In addition, for chaotic systems the action derivatives will 
eventually decouple from the weightings in the correla- 
tion so that in the semiclassical limit the random matrix 
theor y pr ediction is recovered stating that no correlations 
exist [ill [Tgj . On the other hand, for finite h with the 
average computed up to the Heisenberg time, there is a 
possibility of short orbits that do not decouple in chaotic 
systems, so correlations may exist. For near-integrable 
systems there is no decoupling even in the long time limit. 

Combining the previously derived results in chaotic 
systems with the newly derived results for near-integrable 
systems, for mixed phase space systems we attempted 
to separate the orbits for applying the Hannay-Ozorio 
de Almeida sum rule to the different regions of phase 
space. Several criteria were chosen as a basis for the 
separation. All of them led to inconsistencies in the 
interpretation. It seems that all of the unstable orbits 
spend some part of their evolution mimicking stable mo- 
tion and do not give standard diffusion contributions to 
action derivatives. Nevertheless, a different semiclassi- 
cal procedure [24| , not going through the Hannay-Ozorio 
sum rule derivation, gives a theoretical prediction that 
matches well with the level velocity variances indepen- 
dently of the dynamical regime of the system. 



The method 17] to obtain the variance of the level 
velocities begins with the smoothed spectral staircase 



N e (E,X)=J2^( E - E nW) 



(Al) 



Taking the derivative with respect to the parameter, we 
obtain 

3N e (E,X)^ s{E _ Mx)) BEg) (A2) 



dX 



The quantity e is an energy smoothing term which will be 
taken smaller than the mean level spacing. Upon energy 
averaging 



dN e (E,\) 
dX 



= d(E, X) 



dE n (X) 
dX 



(A3) 



where d(E, X) is the mean level density which is the recip- 
rocal of the mean level spacing and proportional to Ti~ d . 
The mean level spacing is related to the phase space vol- 
ume by V = dh d . Thus, in order to obtain information 
about the level velocities, we will evaluate the spectral 
staircase. 

The semiclassical construction of the spectral staircase 
is broken into an average and an oscillating part 



N e (E, X) = N(E, X) + N osc (E, A) 



(A4) 



The average staircase N(E, X) is the Weyl term and to 
leading order in h is given by 



N(E,X) = 



1 

1? 



8(E — H(1, 6; X)dldO 



(A5) 



The oscillating part of the spectral staircase is a sum 
over rational tori with topology M of the unperturbed 
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Hamiltonian, Hq(1) ^| 



N osc (E,X) 



M 



xRe|exp(^ (A) -^ + ^ N 
x[J (A5 M (A)A)-m(A)Ji(A5 M (A)A)] 



x exp 



(A6) 



The phase jyjvi = M rj where r\ are the Maslov indices 
and /?m is equal to the sign of the determinant of the 
curvature matrix. 

The variance of the level velocities is obtained by tak- 
ing the mean square of the counting function derivatives 



dN e (E, A) 



EE 



dE n (X) dE m (X) 



OX J dX d\ 



x S e (E - E n (X))6 e (E - E m {\)) ) 

1 1 



(A7) 



For a nondegenerate spectrum, the summation is nonzero 
only if n = to because of the product of the two delta 
functions. Since Lorentzian smoothing is applied, then 



5l(x) 



2tt£ 



<5 £ /2(^) 



(A8) 



for e d . Thus we have 



dN e (E,X) 
dX 



d f dE n (X) 
2ire \\ dX 



(A9) 



The summation over the rational tori is most sensitive 
to the changing actions and periods because of the asso- 
ciated rapidly oscillating phases, i.e., the division by Ti 
in the exponential. Since the energy smoothing term e is 
taken smaller than a mean level spacing, it scales at least 
by h d and the derivatives of the period vanish as Ti — > 0. 
Thus, only the derivatives of the actions are considered, 
and the oscillating part of the staircase yields 



dN osc (E) 
dX 



fdS 



M 



V dX 



h (d+D/ 



T)7lE' AMlln l- ex P 



M 



[J (s\) - ia(X)Jx(s\)] + i 



ti 

f dAS M 
\ dX 



ia(X) , , 



o.MM-$^)KA10) 



After making the diagonal approximation and energy av- 
eraging, the variance of the level velocities for a near- 



integrable system is 



J 



47re^//,2 J ( dSm \ { , 2( ^,~2 T 2, m , ( 9S M \ ( dAS m 



M 



\ dX 



-(Jo(sx)-J2(sx)) 2 + J?(sx) 



;Jo{sx){Ja{sx) - J2(sx)) + aJ%(s\) 



(All) 



r 



The middle term vanishes since (dSM./dX) and 
(dASm/dX) are uncorrelated and both average to zero. 



APPENDIX B: OVERLAP INTENSITIES 



The variance of the intensities is derived by examining 
the oscillating part of the strength function 



S aosc (E,X) = — Im / (a\0)G(6,6';E)(9'\a)d0de' 

7T 



(Bl) 
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where 



0(0,9'- E) 



/2 



ex-p[iS J (0,0';E)/h 



m/Jtt/2] (B2) 



is the semiclassical energy Green's function. The above 
sum is over all paths that connect to 9' on a given 
energy surface E. D s is a determinate involving second 
derivatives of the actions 



d 2 S 


a 2 s 


aeae' 


aeaE 


a v s 


a r s 


dEdO' 


7)W 



(B3) 



We assume that the Gaussian wave packet in Cartesian 
coordinates can be define locally as a Gaussian in action- 
angle variables using the Wigner function 



A w (I,0) = 2 d cxp{-(AI) 2 /?i 2 - (AO) 2 } 



(B4) 



where (AI)j = <jj(I — l a )j and (A0)j = (9 — O a )j/aj. 
The inverse transform is evaluated by stationary phase 
to obtain 



{a\0)(0'\a) = 



1 



x exp 



(AO) 2 (AO') 2 



(B5) 



i(i-i a )-(e-e') 



The Gaussian localization in action-angle space implies 
that the dominant contributions to the strength function 
arise from tori whose trajectories are closed, i.e., rational 
tori Im- An expansion about the center of the Gaussian 
wave packet yields 



Sj(0, 0'; E) — Syi(O a ,O a \ E) + Ij ■ (0 — 0') 



(B6) 



where quadratic terms have been neglected. Using the 
above expansion of the action, the integrals over the an- 
gles can be performed by stationary phase giving 



^ct,osc 

(E,X) 



2 (d+l)/2-Q a3 

7r i/2 ft («H-D/a 



/2 



(B7) 



M 





r-(Ai) 2 i 


j exp 


h 2 



The above result for the oscillatory part of the strength 
function can be used to obtain the variance of intensities 




exp 



-2(AI) 2 



(B8) 



TABLE I: Fixed Lyapunov exponent cutoff, /i is the fraction 
of orbits that have a Lyapunov exponent less than \i = 0.5. 



k 


C(fc) 


/ 


h 


1.5 


9.62 x 10" 5 


0.982 


1.000 


2.0 


1.11 x 10~ 4 


0.856 


0.798 


2.5 


1.76 x 10~ 4 


0.358 


0.336 


3.0 


3.66 x 10 -4 


0.161 


0.142 


3.5 


3.58 x 10~ 4 


0.166 


0.146 


4.0 


3.67 x 10~ 4 


0.118 


0.107 



the density of states given by Berry and Tabor and the 
density of states is the energy derivative of the spectral 
staircase function, then \D„\ = T^A^/ '(27r) d_1 . Hence, 
given the Hannay-Ozorio sum rule for ^ M A 2 ^, Eq. 10, 



then 



V 



M 



(27r) d 



dT- 



(B9) 



Since the Gaussian weighting of the action variables is 
independent of the period, it follows that 



exp 



-2(AI) 2 



(B10) 



M 



APPENDIX C: BALLISTIC VS. DIFFUSIVE 
BEHAVIOR 

We tried three different methods for separating out 
the ballistic and diffusive orbits. The fraction of ballistic 
orbits / was deduced by 



/ = 



C2(k) 
47T 2 C(fc) 



(CI) 



where 02(h) was determined by Eq. (|44l) . The first 
method was to divide the orbits by a fixed value of the 
Lyapunov exponent. It turned out that a Lyapunov ex- 
ponent of 0.5 gave fairly good result, but there is no 
explanation of this value. Table [I] summarizes these re- 
sults. The next method was to take twice the fraction 
of stable orbits as the number of ballistic orbits. This is 
summarized in Table |n] For the final method we took 
the greatest number of orbits up to some Lyapunov expo- 
nent that did not change £(fc) significantly as the number 
of ballistic orbits and is summarize is Table IIIII 



where we have used the diagonal approximation for the 
sum over M. Since the trace of the Green function is 
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TABLE II: Twice the regular orbits, is twice the fraction of 
regular orbits and fi2 is the cutoff of the Lyapunov exponent 
associated with f'2- 



k 


C(fc) 


/ 


h 


M2 


1.5 


1.17 x 10" 4 


0.812 


0.687 


0.318 


2.0 


1.49 x 10~ 4 


0.638 


0.526 


0.443 


2.5 


2.19 x 10~ 4 


0.288 


0.256 


0.420 


3.0 


2.19 x 10~ 4 


0.269 


0.241 


0.590 


3.5 


2.21 x 10~ 4 


0.268 


0.246 


0.702 


4.0 


2.29 x 10~ 4 


0.189 


0.179 


0.713 



TABLE III: Stable quadratic, fz is the number of orbits in- 
cluded in the ballistic term which does not change its coef- 
ficient significantly. /13 is the corresponding Lyapunov expo- 
nent. 



k 


C(k) 


/ 


h 


M3 


1.5 


1.84 x 10 -4 


0.516 


0.412 


0.2 


2.0 


2.58 x 10~ 4 


0.368 


0.302 


0.3 


2.5 


3.30 x 10~ 4 


0.191 


0.165 


0.2 


3.0 


4.26 x 10~ 4 


0.139 


0.124 


0.4 


3.5 


4.01 x 10~ 4 


0.148 


0.137 


0.4 


4.0 


4.31 x 10~ 4 


0.105 


0.100 


0.4 
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